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ABSTRACT 



The evidence for a rotation of the e Eridani debris disc is examined. Data at 
850 (a.m wavelength were previously obtained using the Submillimetre Common User 
^ ' Bolometer Array (SCUBA) over periods in 1997-1998 and 2000-2002. By fitting 

after shift and rotation operations, images from these two epochs were compared to 
recover proper motion and orbital motion of the disc. The same procedures were then 
^ . performed on simulated images to estimate the accuracy of the results. 

' Minima in the plots indicate a motion of the disc of approximately 0.6" per year 

00 \ in the direction of the star's proper motion. This underestimates the true value of 1" 

lO ' per year, implying that some of the structure in the disc region is not associated with e 

' Eridani, originating instead from background galaxies. From the fitting for orbital 

, motion, a counterclockwise rotation rate of ^2.75° per year is deduced. Comparisons 

with simulated data in which the disc is not rotating show that noise and background 
galaxies result in approximately Gaussian fiuctuations with a standard deviation ±1.5° 
r"| . per year. Thus counterclockwise rotation of disc features is supported at approximately 

i' a 2-(T level, after a 4-year time difference. This rate is faster than the Keplerian rate 

' \ of 0.65° per year for features at ?a65 AU from the star, suggesting their motion is 

^ . tracking a planet inside the dust ring. 

' Future observations with SCUBA-2 can rule out no rotation of the e Eridani dust 

, clumps with ~4(T confidence. Assuming a rate of about 2.75° per year, the rotation of 

the features after a 10-year period could be shown to be ^1° per year at the 3cr level. 



Key words: Kuiper Belt - circumstellar matter - planetary systems: formation 
planetary systems: protoplanetary discs - submillimetre. 



1 INTRODUCTION 

1.1 e Eridani: A Special Case 

Debris discs around nearby stars represent extra-solar ana- 
logues of the Kuiper Belt. It is in this context that study- 
ing the e Eridani dust ring is of particular interest since 
it is of spectral type K2V, not too dissi milar to the Sun , 
but with a mu ch you nger age of 0.8 Gyr llSong et al.ll2000l: 
iDi Folco et"al] l2004ll and thus represents an analogue to 
the young solar system. It is generally thought that proto- 
planetary discs evolve into debris discs after planets have 
formed, a proces s that is completed within 10-100 My r 
after star birth iSchutz et alj l2004l : iHoUand et all 11993) . 
iGreaves et alJ i2004l) have noted that there is an apparent 
lack of overlap between stars with debris discs detectable 
at submillimetre wavelengths and those with radial veloc- 
ity planet detections. There are only a few stars with pos- 
itive detections for both an infrared dust-excess and a ra- 
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dial velocity planet (iDominik et alj|l998l : iBeichman et alj 
I2OO5I) . While this is likely to reflect the small number 
of stars searched for both phenomena, e Eridani remains 
unique in having resolved structure within the debris ring 
iGreaves et alj Il998l . l2005l) pl us an inner gas gia nt planet 
detected by Doppler wobble jHatzes et al.l l200(]|) . Further 
evidence for this pl anet comes from the forced offset of the 
centre of the ring llGreaves et al.ll2005l) . consistent with a 
forced eccentricity of the dust particles (.Wvatt ct al..i.l999i) . 
The ring structure also appears perturbed as if by a more 
distant gas giant, but this body has not yet been imaged, 
with a maximum mass constr aint at around the dus t ring 
radius of approximately 5 Mj iMacintosh et al]l2003l) . 

1.2 Planet Hunting by Tracking of Disc Features 

An infrared excess around e Eridani was first detected dur- 
ing photometric measurements by the Infrared Astronomical 
Satellite (IRAS) llAumannlll9'83l . soon after the initial dis- 
covery of an IR excess around Vega iAumann et alJfl984^ . 
Subsequent submillimetre observations at 850 \ixn with the 
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Su bmillimetre Common User Bolometer Array (SCUBA) 
by iGreaves et alj il998l) resolved the disc and showed a 
ring-like structure. The faint ring surrounding e Eridani 
has not yet been suc cessfully imaged at optical wavelengths 
JProffitt et alj|2004) . 

Further SCUB A observations made in 2000-2002 
JCreaves et alJl2005^ allow us to study the motions of the 
substructure in the ring over a 5 year period. If clump fea- 
tures are tracking the orbital motion of a planet within the 
ring, they should appear to rotate faster than the Keplerian 
rate at the ring radius; this would be an unambiguous signa- 
ture of a planet with the forced period providing a measure 
of its orbital semi-major axis. Clumps will appear in a disc 
where planet formation has taken place because the plan- 
ets will migrate outwards via angular momentum exchange 
with the remaining planetesimals. As the planets migrate 
outwards, the associated gravitational resonances also move 
out, thus trapping dust parti cles outside the planet's orbit 
into mean motion resonances (IWvatt et al.ll999l . l2003f) . The 
Kuiper Belt in the Solar System probably formed in a simi- 
lar way, where the material is thought to have been pushed 
outwards by a mean motion resonance associat ed with Nep- 
tune as it migrated further away from the S un iGomesl2003l : 
iLevison fc MorbideUill20o3 : lMalhotrallT995l) . 

Circumstellar dust, trapp ed into resonances by an 
outwardly m igrating planet iWvatt et alJ Il999l. 120031: 
lOuillen fc Thor ndikc 2002). is expected to be observable at 
sub-millimetre wavelengths. Whilst the Vega and Fomalhaut 
disc asymmetries ha ve been modelled llWvatt et aljl20o3 
IWvatt fc Dentl2002^ ■ the e Eridani dust ring offers a unique 
prospect to directly track motion of clumps over time, be- 
cause of its favourable inclination (~25° from face-on to the 
observer) and the fact that this is one of the closest stars 
to the Sun at a distance of 3.2 2 pc. Models o f the distri- 
bution of the e Erid ani clumps llOuillen fc Thor ndikc 2003 
lOzernov et"alll200(f) suggest the perturbing planet lies «40- 
60 AU from the star, thus the clumps are expected to orbit 
with a period of ~280-520 years. These rotation rates of 
0.7-1. 3° per year, c orresponding t o particular resonan t pe- 
riods jQuillen fc Th orndikc 2002: lOzernov et aljlioooll . re- 
quire long times to detect. Here we present a preliminary 
analysis over 5 years of data, making use of image fitting to 
recover statistically any small rotation present. 



2 OBSERVATIONS 

Th e results used here co nsist of the 5 year dataset presented 
bv ICreaves et all feOOST) . Observations of e Eridani at 850 
^.m were made with SCUBA between 1997 August and 2002 
December, in total comprising 56 images and an integration 
time of 33.5 hours. The dataset is split into two parts, one 
including data taken during 1997-1998 and the other from 
2000-2002; the observing runs were bunched in time so that 
the effective mid-points of the two periods are around 4 years 
apart. The 2000-2002 dataset has a noise level lower by a 
factor of two, as a result of longer duration and better sensi- 
tivity. Data were also obtained at 450 vim, but the SCUBA 
filter used in 1997-1998 had a low throughput, and so this 
image can not be used for time-dependent studies. 

The analysis made here is of the co-added 1997-1998 
data compared to the co-added 2000-2002 results. The stel- 



lar proper motion is -1" in RA per year, so any disc emission 
associated with the star should shift by approximately 4" 
west, with some blurring because the observing periods were 
spread out. The diffraction-limited FWHM beam size at 850 
\im was 15" by 15.5", but the data have bee n smoothed us- 
ing a 7" Gaussian to an effective 17" beam dCreaves et alJ 
120051) . The net pointing errors are expected to be small, be- 
low the level of the 1" cell size in the images. The noise for 
the 1997-1998 dataset is 4.3x10"^ mjy arcsec"^ (0.96 mjy 
beam~^) with 2.4x10"^ mJy arcsec"^ (0.54 mJy beam^^) 
for the 2000-2002 dataset. Photospheric emission of 1.7±0.2 
mJy has been subtracted from the images, for ease of com- 
parison with simulations that neglect the star. 



3 SIMULATED DATA 

The proper motion and rotation of dust features associated 
with the star can be assessed by translating and rotating the 
1997-1998 and 2000-2002 co-added images to find the best 
match. The accuracy of the results was evaluated by compar- 
ison with simulated data with comparable noise, background 
galaxies and a foreground consisting of clumps embedded in 
a smooth ring of radius 20" centred on the star. 

Each simulated image was constructed from a set of 
frames, each representing an individual observation taken 
with SCUBA. The simulated analogue of the 1997-1998 
dataset comprised 22 frames of equal depth, which when co- 
added reproduced the noise in the actual data at the final 
17" resolution. The range of proper motion offset compared 
to the last real data taken was 4-5". A similar procedure was 
followed to simulate the 2000-2002 data, using 34 frames and 
a proper motion range of 0-2.5". The difference in depths 
between the real 1997-1998 and 2000-2002 datasets was ac- 
counted for by both the number of images in each simulated 
dataset and the noise levels that were input into each set. 
While the number of input frames is the real value, the equal 
noise value per frame is a simplification (as it was too com- 
plex to simulate the duration and observing conditions of 
every real frame). Thus the effective mid-points of the two 
observing periods, about 4 years apart, may not be exactly 
matched in the simulated results. 

Each simulated frame representing an observation by 
SCUBA was constructed from random noise, random back- 
ground galaxies and a foreground (ring and clumps embed- 
ded in the ring) of comparable brightness to the observed 
data. The same background galaxies were used in both the 
1997-1998 and 2000-2002 frames but a different sample gen- 
erated for each simulation. The chopping motion of the sec- 
ondary mirror of the telescope used to determine sky levels 
was also simulated. All of these effects contribute to the real 
images; the data should thus represent one possible outcome 
of the simulations, to within the accuracy of the simplifica- 
tions used. 



3.1 Simulation details 

The noise image was made by choosing the fiux for each in- 
dividual pixel randomly from a Gaussian distribution and 
then smoothing spatially with a 7" Gaussian as was done 
for the observed data. This assumes the noise per 1" pixel 
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is statistically independent. The mean for the noise distri- 
bution was zero and the a value was chosen so that the final 
simulated images had a standard deviation of flux per beam 
matching the real dataset once the galaxy population had 
been added. 

The integral galaxy counts N(>S) for each flux, S, at 
850 \ ira were modell e d with Poisson statistics using numbers 
from lBarnard et al.l i2004l) . by a double power law given by 

iV(> S) oc S"" (1) 

with a=0.94 for S < 1.18 mjy and a=2 for S > 1.18 mjy. 
Each galaxy's point-like flux was smoothed with a 2-D Gaus- 
sian to reflect the size of the beam of the JCMT at 850 |J.m 
and then placed at random coordinates chosen from a uni- 
form distribution. Thus no area of the image was favoured 
and any possible galaxy clustering was neglected. 

The foreground was modelled by adding beam-sized 2- 
D Gaussian regions of flux to a ring. The ring was created 
by placing Gaussians at every pixel at distances between 
17.5"-22.5" from the star, producing an ann ulus centred on 
the s tar with the radius and width observed (iGreaves et alJ 
I2OO5I) . The flux per pixel in the ring was chosen to match the 
smooth l evel observed ( i .e. be tween clumps) of f»10^^ mJy 
arcsec~'^. lGreaves et alj (|200^ identifled three clumps with 
possible rotation located northeast, northwest and southeast 
of the star at a radius of 20" and thus three clumps were 
added at these positions in the simulated ring. The peak 
fluxes of the clumps were set at f»1.5 x lO^'^ mJy arcsec^^ 
so that the total foreground had peak fluxes at the clump 
locations of «2.5 x 10~^ nijy arcsec"^. This is the mean of 
the total fluxes towar ds the three candida te moving clumps 
in the observed data (iGreaves et aLllioOsh . 

Proper motion of the disc was simulated by translat- 
ing the foreground relative to the background for each of 
the simulated observations, by a distance corresponding to 
the epochs of individual real frames. In the observed data, 
a correction for the annual proper m otion of e E ridani {^a 
= -0.976" yr-\ /i^ = 0.018" yr'^) JESA||1997I^ was made 
by sorting an d shifting the frames to a precision of 0.5" bins 
in R.A. only dGreaves et al.ll2005^ . This correction was also 
performed on the simulated observations, and so the back- 
ground galaxies (fixed with respect to the sky) will appear to 
move at the rate of the proper motion but towards positive 
R.A., in a co-ordinate frame co-moving with the star. 

Also simulated was the chopping procedure performed 
by the JCMT, used to remove sky fluctuations. 'Blank' re- 
gions of sky on either side of the held of interest are observed 
interleaved with the on-source observations, and subtracted 
to leave onl y astronomical signal; f or comments on the lim- 
itations see lArchibald et al.N2002) . If background galaxies 
lie within the SCUBA field of view at either off-source posi- 
tion, "holes" with magnitude equal to half the brightness of 
the galaxies being chopped on to will appear in the final ob- 
served image. Due to variations in the chop direction, such 
holes will be blurred in the on-source residual frame and thus 
difficult to recognise. As chopping introduces extra fluctua- 
tions and modifies the background galaxy contributions, it 
was important to include it in the simulations. 

The noise, background galaxies and foreground images 
were constructed for a 400 by 400 pixel square frame but 
the field of view of SCUBA on the JCMT is roughly circular 
with a radius of 72" . It was necessary to simulate this larger 



area of sky to allow the chopping procedure to be simulated. 
Flux values for each pixel within an area of sky were read 
into 2-D arrays at the position of the source, M, and the 
two off-source positions to the left, L, and right, R, of the 
source. The pixel fluxes used in the flnal images, I, were then 
calculated using 

/(.,y)^M(x,,)- (^(-'^);^("'^» . (2) 

Chopping in azimuth was used for the real data, so that 
the left and right chop throws fell on different sky regions 
depending on the source elevation. This was accounted for 
in the simulations by varying the angle between the chop 
throw and the line of constant right ascension (horizontal in 
the images) randomly between -50° and 50°. In the obser- 
vations, chop throws of 80" (5% of images), 100" (55%), and 
120" (40%) were used. A chop throw of 120" was used for aU 
of the simulated chopping. 

Some small effects were neglected in the simulations. 
The disc is fairly close to face-on to the observer, but in- 
clined at ~25° from the sky plane and so actually appears 
as an ellipse in the sky. Assuming a rotation rate of 2.75° 
per year, a clump moving along the direction of the minor 
axis of the ellipse will have a motion that is only ^^0.3" less 
than a clump moving parallel to the direction of the major 
axis. Therefore, it was reasonable to approximate the ring 
as circular to simplify the rotation of the foreground image. 
The annual shifts in position due to parallax for e Eridani 
are ±0.3" and these were neither accounted for in the tele- 
scope tracking software nor included in the modelling. In 
the observed data 68% of the images co-added were taken 
between August-October and so the effect of smearing due 
to parallax errors will be less than 0.3". The dust ring was 
also simulated as being cent red on the star , negl ecting the 
offset of 1.5-2" identified bv IGreaves et all i2005ft and sus- 
pected to be forced by the inner planet. We assume the bulk 
motions of the clumps still follow circular orbits. 

3.2 Simulation output 

The real images are shown in Figure 1 and examples of the 
simulated images are shown in Figure 2. The simulated im- 
ages recreate a ring with comparable brightness, size and 
morphology to that seen in the observations. The regions 
surrounding the ring also show comparable levels of noise 
and density of background objects. Differences between the 
two simulated images are at a similar level to those in the 
real data of Figure 1. Although the image structure looks 
complex, this is in part because the false-colour scale both 
includes a wide fiux range (negative as well as positive sig- 
nals) and also tends to enhance features with low levels of 
contrast, as discussed bv IGreaves et all (|200^. 

The simulation procedure was also repeated to resem- 
ble a future observation that could be taken in 2007 with 
SCUBA-2 - a replacement camera for SCUBA with greater 
sensitivity, fidelity and field of view llAudlev et all l2004l: 
iHolland et ai]|2003^ . For this simulation, the noise was set 
to zero as it will fall below the background confusion in 
an observation of only ~1 hour. The chopping part of the 
simulated data algorithms was also disabled since SCUBA- 
2 will sample rapidly enough that sky fluctuations will be 
negligible and so chopping will not be required. (Chopping 
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Figure 1. Observations of the e Eridani debris disc taken witli SCUBA on tlie JCMT. (Left) 1997-1998 dataset, (averaged over 22 images) 



in a 1" pixel grid in R.A., Dec coordinates (north is up, east is left), plotted with the dynamic ranges min= 
mjy arcsec"^. (Right) 2000-2002 dataset (averaged over 34 images) with the same parameters. 



-1x10" 



=2.5x10- 





Figure 2. As for Figure 1 but for simulated SCUBA images with no rotation of the foreground clumps: (left) example of simulated 
1997-1998 image, (right) the corresponding 2000-2002 image. Other parameters arc as in Figure 1. 



can be re-introduced in the data analysis where necessary 
for matching SCUBA-2 images of e Eridani with those al- 
ready taken with SCUBA.) Any candidate rotation of the 
disc identified here can thus be checked by extending the to- 
tal timeline of observations to 10 years (from 1997 to 2007) . 



4 DETECTING A ROTATION 

The simplest method of measuring rotation would be to 
track individually i dentified clumps. However, clump ex- 
tract ion algorithms dWilliams et alj[l994l : lBertin fc Arnoutj 
Il996l) are not well adapted to finding clumps and arcs within 



a ring. Hence a x^-fitting technique was adopted instead, 
with the advantage of using all the information in the im- 
age simultaneously. The per pixel value for two images, 
Ii and I2, with measured noise a\ and 02, using only pixels 
in a box with its lower left and upper right corners defined 
by (x,y) cartesian co-ordinates (ai,fei) and (02,62), can be 
calculated using 

2 ^ {h{x,y)-h{x,y)f 
^ (a2-ai + l)(&2-6i + l)K + a2)- 

The per pixel values between the 2000-2002 dataset 
and the 1997-1998 dataset were calculated as the 1997-1998 
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Figure 3. P^r pixel curves for observed data and simulations. (Left) Fits for annual proper motion shift from the 1997-1998 to the 
2000-2002 dataset. Thick curve shows observed data; other curves are 5 random simulations with no proper motion or rotation. Positive 
shifts are to the right (-ve in R.A.) and the arrow indicates the known proper motion of e Eridani = -0.976" yr~^). Small deviations 
in the curves are caused by slight imperfections in the linear interpolation. The narrow curve (Sim D) originates from a very bright galaxy 
in the ring, forcing a minimum that must be close to zero. (Right) Fits for annual rotation rates, comparing 1997-1998 and 2000-2002 
datasets and with corrections made for proper motion. Thick curve shows observed data; thin curves are 5 random simulation pairs with 
proper motion corrections included but no rotation. 



dataset was shifted in right ascension and rotated relative to 
the 2000-2002 dataset, using the tasks "SL IDE" and "RO- 
TAT E" in the KAPPA software package dCurrie fc Berrvl 
|2004) . The minima in plots of x'^ per pixel as a function 
of right ascension shift or rotation angle then identified the 
best solutions for each (independently). The region of the 
image included in the x^ fit was defined by a 70"-square 
box, as this included all of the disc with as little off-source 
background as possible; extra background galaxies included 
in a larger box tend to cause the x^ fits to incorrectly favour 
a shift and rotation of zero. 

For the simulated SCUBA-2 observation, the x'^ per 
pixel values were calulated twice, once for the future 2007 
observation compared to the 1997-1998 SCUBA dataset, and 
again for 2007 compared to the 2000-2002 SCUBA data. 
The two x'^ values corresponding to the same shift or ro- 
tation rate were then averaged to create the final x^ per 
pixel curves. We also performed simulations with no galax- 
ies or noise, and verified that the true rotation values were 
recovered correctly. 

There are three possible ways in which the regions of 
emission in the ring may be behaving: 

(i) The ring consists of only background galaxies and 
therefore the features are not tracking with the proper mo- 
tion of the star or revolving around the star's position. 

(ii) The ring features are tracking with the proper motion 
of the star but are revolving around the star at the Keple- 
rian period of the ring which would not require a planet to 
explain the motion. 

(iii) The ring features are tracking with the proper motion 
of the star and revolving around the star at the period of a 
planet at tens of AU. 

To establish which of these explanations gives the best 
description of the data, the observed data were fitted for the 
proper motion of emission in the disc region, and then fitted 
to establish a rotation rate of the ring features. Then for 



comparison, 100 sets of simulated data were made with no 
proper motion or rotation of the foreground, to quantify the 
random fluctuations due to noise, background galaxies and 
chopping on to background galaxies. This was then repeated 
but assuming the foreground was tracking with the proper 
motion of the star and with the foreground rotated at rates 
of 0° (i.e. no rotation), 1°, 2°, 4° and 10° per year (unreal- 
istically large but included for comparison). The purpose of 
trying different rotation rates was to examine how well the 
true rate is recovered, given that non-moving background 
galaxies introduce a bias towards a x^ minimum around zero 
degrees. 



5 RESULTS AND DISCUSSION 

Figure 3 shows x'^ per pixel as a function of the annual 
proper motion and rotation rate of the 1997-1998 image rel- 
ative to the 2000-2002 data. 5 examples of simulations are 
shown for comparison. Table 1 lists the statistics from the 
full sets of simulations, each set comprising 100 comparisons 
of pairs of images. 

5.1 Motion of the Disc 

Figure 3 (left panel) shows that in the observed data, the 
emission in the ring region has shifted at a rate of 0.6" per 
year to the right (negative in R.A. ) when tracked from the 
1997-1998 to the 2000-2002 epoch. lOreaves et alJ lj200.'jl ar- 
gued that at least some of the disc features are tracking 
with the proper motion of the star, i.e. 1.0"/yr to the right. 
The fact that the net motion measured is smaller may be 
explained by the presence of stationary background features 
behind the ring, which by themselves would produce a min- 
imum in the x'^ curve at 0" per year. A result intermedi- 
ate between 0"/yr and l"/yr thus indicates a mixture of 
tracking and stationary features. The net value of «0.6"/yr 
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Table 1. Proper motion shifts and rotations from fitting of 6 simulated SCUBA datasets each consisting 
of 100 simulated 1997-1998 and 2000-2002 images with foreground background galaxies, noise and chopping 
included. 



Simulation 


Simulated 


Average 


Motion 


Simulated 


Average 


R.otation 


So- Clipping 


Average 


Set No. 


Motion'^ 


Motion 




Rotation 


Rotation 


cr 


Results 


Min Value 




(" yr-i) 


{" yr-^) 


{" yr-i) 


(° yr-') 


{° yr-^) 


r yr-^) 


(No. Excluded) 




1 





0.01 


±0.28 





-0.03 


±1.23 


5 


1.92 


2 


+1 


0.79 


±0.38 





-0.15 


±1.58 


4 


2.03 


3 


+1 


0.78 


±0.42 


-1 


-0.68 


±1.83 


5 


2.03 


4 


+1 


0.73 


±0.44 


-2 


-0.96 


±1.87 


2 


2.06 


5 


+1 


0.65 


±0.47 


-4 


-1.92 


±2.55 


76 


2.14 


6 


+1 


0.63 


±0.47 


-10 


-4.83 


±13.71 


_c 


2.42 



Positive values signify a direction negative in R.A. 
^ 4cr clipping used for these values 

a clipping not possible due to broad distribution of results 



suggests that the majority of features in the ring region are 
truly associated with e Eridani. 

The simulations (Table 1) reproduce this result closely. 
An input proper motion of l"/yr is recovered as on aver- 
age «0.7"/yr of net motion (with little dependence on any 
assumed ring rotation). Thus the stationary galaxies do in 
fact reduce the net motion measured in the same sense as in- 
ferred in the real data. The standard deviation of the proper 
motion solutions is approximately ±0.4"/yi'- Thus the real- 
data solution of 0.6"/yr, differing by only ~0.l"/yr from the 
typical simulation result, is well within the scatter. 

Further tests confirmed the validity of the simulations. 
When the observed data were corrected for the proper mo- 
tion of the star, the solution obtained was 0.4" per year to 
the left (positive in R.A.). This is a measure of the amount 
of emission that is fixed in the sky and so moving to the 
left in the co-ordinate frame of the star. This is the value 
expected, as the difference between the l"/yr if all the emis- 
sion moved with the star and the 0.6" actually measured. 
A set of simulations was also performed to check that a 
non-moving disc returned the correct result of 0" per year. 
The measured mean was very close to zero (Table 1) with 
a standard deviation of ±0.25"/yr. This is a measure of the 
fluctuations that can occur due to random noise and chop- 
ping onto background galaxies; it is somewhat smaller than 
the ~0.4" obtained in the more complex cases with both sta- 
tionary and moving flux contributions. Finally, we verified 
that the small 'ringing' effect seen in the curves in Figure 3 
does not affect the position of the minimum in the ^ curves. 
This ringing arises from a small jump at integer values, ap- 
parently in the bilinear interpolation of the "SLIDE" task 
in KAPPA. Plotting the curves using only integer values 
changed the average motions recovered by only '-^O.Ol". 

5.2 Rotation of the Disc 

The fit between the observed 1997-1998 and 2000-2002 
datasets (corrected for proper motion) shows a rotation of 
«2.75° per year counterclockwise (see Figure 3, right panel). 
This is in the same direction as but somewhat larger tha n 
the rotation of '-^1.5° suggested bv iGreaves et alT i2005l) . 



It is also lar ger than the rotation rate of 0.7°-1.3° per year 
predic ted bv lOzernov et all i2000l) and lQuillen fc Thorndikj 
l|2002), assuming dust is trapped in resonances with a planet 
at f»40-60 AU from the star. 

The simulated data are used as a guide to what rota- 
tion rates can be confirmed or ruled out. In the simulation 
results (Table 1), the x^ fit results from sets of 100 image- 
pairs showed that the rotations obtained were approximately 
Gaussian distributed, for the cases where the rotation rate 
was 0°, 1° and 2° per year. There were, however, a few spu- 
rious results which involved an unphysical disc feature rota- 
tion rate (total rotations of >50°). These results occurred 
when a bright clump was randomly generated by a com- 
bination of noise and background galajcies and fortuitously 
positioned in the dust ring to give a deeper minimum in the 
X^ curve corresponding to a large rotation angle. These out- 
comes (^ 5 per set of 100 simulations) did not fit into the 
Gaussian distribution and were removed by clipping at the 
5(7 level. After removing these points, the sample standard 
deviation was typically ~1.5° of rotation per year. For ro- 
tation rates of 4° and 10° per year, the rotations obtained 
followed a broader distribution, and so no reliable measure 
of the rotation was obtained. 

The true rotation rate is generally not recovered in the 
simulations, with an average annual rotation measured gen- 
erally about half of the input rotation. This can be under- 
stood as the effect of the stationary galaxies, which will in- 
troduce a bias towards a solution of zero. (When no rotation 
is input, a solution close to zero is in fact recovered.) Only 
simulations in roughly the -la tails of the distributions re- 
cover the true input rotation rate (Table 1). There is there- 
fore a moderate probability that the rate measured in the 
real data, of «2.75° per year, is actually an underestimate 
for the true value. 

Table 1 shows that for the simulated data with no ro- 
tation of the foreground, the sample standard deviation was 
1.44° per year. As the measured rotation was 2.75° per year, 
this means that the null hypothesis (zero or very slow ring 
rotation) is ruled out at the w2cr level. However, the actual 
rotation rate can not be extracted from comparison with the 
simulations as the results overlap; that is, the standard de- 
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Table 2. Proper motion shifts and rotations from fitting of 6 simulated datasets including the SCUBA 
frames of Table 1 and a SCUBA-2 image to be obtained in 2007 (for which no noise is included). 



Simulation 


Simulated 


Average 


Motion 


Simulated 


Average 


Rotation 


So- Clipping 


Average 


Set No. 


Motion*^ 


Motion 


(7 


Rotation 


Rotation 




Results 


Min Value 




(" yr-i) 


{" yr-i) 


{" yr-i) 


(° yr-i) 


(° yr-i) 


(° yr-i) 


(No. Excluded) 




1 





0.00 


±0.07 





-0.05 


±0.26 





2.05 


2 


+1 


0.81 


±0.21 





0.05 


±0.64 





3.07 


3 


+1 


0.78 


±0.22 


-1 


-0.52 


±0.68 





3.19 


4 


+1 


0.75 


±0.22 


-2 


-1.04 


±1.08 


1 


3.11 


5 


+1 


0.76 


±0.25 


-4 


-2.48 


±1.79 


1 


3.44 


6 


+1 


0.76 


±0.25 


-10 


-3.70 


±6.00 


b 


4.52 



Positive values signify a direction negative in R,.A. 
^ (7 clipping not possible due to broad distribution of results 



viations are wider than the changes in mean value between 
sets with different input rotation rates. Some estimate can 
be made of how plausible a measured rotation rate of 2.75° 
anticlockwise is within each simulation set. We find that the 
probability of measuring at least this large a counterclock- 
wise rotation is 10-25% for rotation rates of 1-4° per year. 

Tests were made to see if any further observational con- 
straints reduced the simulation parameter space. In the real 
data, the brightest signal is 10 mjy within one beam. Galax- 
ies brighter than this within the 70" box are predicted to oc- 
cur in only ~5% of the simulations. Eliminating these might 
reduce the tendency towards solutions of zero rotation, but 
in fact an increase of approximately 20% was the largest 
change to any recovered rotation rate when this was done. 

We also verified that small defects seen in the curves 
about zero probably arising from interpolation in the rota- 
tion process did not significantly change the results. Removal 
of the defects only changed the rotation rate recovered of the 
non moving non rotating case and only by ~0.001° per year 
and the associated standard deviation by only ^-^0.004° per 
year. 



5.3 Future Observations 

Table 2 shows the results for the simulated datasets with 
a third image with no noise included. This is intended to 
refiect observations that could be made with SCUBA-2 at 
850^tm in 2007, with noise that is negligible compared to ran- 
dom background galaxies. Chopping was re-introduced to 
allow matching of SCUBA and SCUBA-2 images, although 
testing showed that this only sigificantly reduced the random 
fiuctuations in the minima in the 4° and 10° per year cases. 
The minima were found by comparing this simulated 
SCUBA-2 dataset with the two simulated SCUBA prede- 
cessors, as described above. The projected 10-year timeline 
should refine the solutions for the disc motions, and so the 
annual errors are expected to be reduced. 

The proper motion is now recovered at similar mean 
values to before {~0.7" per year) but with lower standard 
deviations. The errors of « ±0.25" /yr are similar to the 
SCUBA-only non-moving simulation, suggesting this is a 
limit set by the fluctuations within the two SCUBA datasets 



from 1997-1998 and 2000-2002. The input rotation is not re- 
covered more reliably than before but the dispersion of the 
results is smaller, significantly so in the cases of up to 2°/yr 
rotation. The number of doubtful results with very large ro- 
tations is also reduced (Table 2). However, the distributions 
of rotation values still overlap between different input rates, 
and so addition of a third-epoch SCUBA-2 image would not 
uniquely identify the true rotation rate. In fact, eliminat- 
ing noise by comparing two hypothetical future SCUBA-2 
observations, with mid-points 4 years apart, still did not re- 
solve this ambiguity. We conclude that fitting alone can 
not identify an orbital period; the paradox inherent in try- 
ing to fit both moving and non-moving components with one 
rotation rate is likely to be responsible. 

With the inclusion of the SCUBA-2 epoch, the case of 
no rotation could be ruled out at the «4(t level. The aver- 
age rotation rate recovered was +0.04° per year with a sam- 
ple standard deviation ±0.74°. Hence if the rotation rate of 
—2.75° per year estimated from the present data were to 
persist, this would be around the -4cr bound of this simula- 
tion. Further, a rotation rate of ^l°/yr would be confirmed 
at around the 3cr level: that is, a continuing measurement 
of —2.75° per year would be just within the -Scr tail of the 
simulation with an input rate of l°/yr. This would support 
clumps orbiting faster than the Keplerian rate at «65 AU 
radius, which is 0.65° per year for a 0.9 Mq star. 

It may also be possible to observe at 450 \mi with 
SCUBA-2 and the improved resolution of 8" at this wave- 
length may allow the rotation to be constrained more accu- 
rately. This would require two epochs of SCUBA-2 data a 
few years apart, and also a better knowledge of the galaxy 
counts at 450 \im. 



6 CONCLUSIONS 

We have identified a motion of the e Eridani disc consis- 
tent with features tracking in the same direction as the 
proper motion of the star. The best fitting proper motion p 
= 0.6"±0.4" per year is most probably less than the star's 
proper motion of 1" per year, confirming that some of the 
features in the vicinity of the disc are likely to be background 
galaxies. The measured rotation rate of 2.75° per year coun- 
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terclockwise is in the same direction but ne arly twice as large 
as that suggested bv lGreaves et alj (|200^; it could also be 
an underestimate as stationary galaxies tend to bias the ro- 
tation solution below the real amplitude. Comparisons with 
simulated data show that the measured rotation is signifi- 
cant at the ~2(T level, and a future image with SCUBA-2 in 
2007 could rule out no or very slow rotation at the ~4cr level. 
The technique of fitting is inherently limited by trying to 
find one solution that fits both stationary galaxies and mov- 
ing ring features, hence a method that identifies individual 
clumps or a method capable of using the proper motion to 
distinguish the foreground and background images will be 
needed to measure the true rotation rate. 

Radial velocity detections of a planet are restricted to 
finding planets out to only a few AU from the star. In the 
case of £ Eridani, finding a planet out to tens of AU via ra- 
dial velocity measurements will be extremely difficult since 
the period of the planet is deduced to be >100 years, and the 
magnitude of the Doppler wobble is reduced by the nearly 
face-on inclination (~25°) of the disc. However, this view- 
ing angle also means that e Eridani provides a unique op- 
portunity to track the motion of the substructure within 
the disc. This is the first ever analysis attempting to track 
clumps rotating in a dusty disc. Future observations with 
SCUBA-2 can confirm the rotation of clumps in resonance 
with a planet at tens of AU, while imaging with submillime- 
ter interferometers such as SMA and ALMA could radically 
shorten the timescale to identify such motions. 
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